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Abstract 



We present a measurement of 7r~/ir + , K~/K + and p/p from p + p collisions at 
y/s = 200 GeV over the rapidity range < y < 3.4. For px < 2.0GeV/c we see 
no significant transverse momentum dependence of the ratios. All three ratios are 
independent of rapidity for y < 1.5 and then steadily decline from y ~ 1.5 to 
y ~ 3. The ir~ /ir + ratio is below unity for y > 2.0. The p/p ratio is very similar 
for p + p and 20% central Au+Au collisions at all rapidities. In the fragmentation 
region the three ratios seem to be independent of beam energy when viewed from 
the rest frame of one of the protons. Theoretical models based on quark-diquark 
breaking mechanisms overestimate the p/p ratio up to y < 3. Including additional 
mechanisms for baryon number transport such as baryon junctions leads to a better 
description of the data. 

Key words: proton collisions, particle ratios, forward rapidity, limiting 
fragmentation, baryon junctions 
PACS: 25.75.q, 25.40.-h, 13.75.-n 



1 Introduction 

The ratios of particle production in hadronic interactions are important in- 
dicators of the collision dynamics [1]. By comparing large and small systems 
over a wide range of phase space, one can address both reaction mechanisms 
in simpler systems and the properties of hot and dense nuclear matter in 
large systems. A thorough understanding of p + p collisions at ultrarelativistic 
energies is necessary both as input to detailed theoretical models of strong 
interactions, and as a baseline for understanding the more complex nucleus- 
nucleus collisions at RHIC energies. Soft particle production from ultrarela- 
tivistic p + p collisions is also sensitive to the flavor distribution within the 
proton, quark hadronization and baryon number transport. Extensive data 
exist near midrapidity, but less is known about the forward rapidity region 
where fragmentation and isospin effects are important. 

In this Letter we present measurements of like-particle charged hadron ratios 
from p+p collisions at a center-of-mass energy of \/s = 200 GeV as a function 
of rapidity y = 0.5 • ln((E + p z ) j[E — p z )) and transverse momentum pt, and 
make a comparison with similar BRAHMS results from the 20% most central 
Au+Au collisions at the same energy. We show that the p + p and Au+Au 
results on pion, kaon and proton like-particle ratios are consistent over three 
units of rapidity, in spite of the expected large differences in dynamics between 
these systems. 

In p + p collisions at RHIC energies two main mechanisms for particle pro- 
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duction are expected. At midrapidity the Bjorken picture [2] predicts that 
particles will be formed mainly from string fragmentation, yielding values of 
antiparticle-to-particle ratios close to unity. At forward rapidities, close to the 
beam rapidity (y b = 5.3 at y/s = 200 GeV), cross-sections are instead known 
to be dominated by leading particles and projectile fragments (the fragmen- 
tation region). This means that the conservation of charge and isospin will 
become increasingly important for particle production as one approaches ?/&. 
The present data on it~ /it + , K~ / K + and p/p show that in p + p collisions at 
y/s = 200 GeV there is a midrapidity region extending out to y ~ 1.5 where 
the particle ratios agree with the Bjorken picture. Above this point the ratios 
start to decrease, indicating the onset of fragmentation region physics. Shifting 
the ratios by the beam rapidity and comparing to lower energy data, we find a 
broad rapidity range where ratios of like-particle production are independent 
of the incident beam energy when viewed from the rest frame of one of the 
protons (limiting fragmentation [3]). 

The traditional quark-diquark breaking picture of a p + p collision fails to 
reproduce baryon transport in available midrapidity data, which has been 
taken as evidence for several additional mechanisms being important at higher 
energies [4,5,6,7]. In this Letter we provide a comparison of different model 
predictions with experimental data, which, especially away from midrapidity, 
provides new constraints for calculations. We show that the commonly used 
event generator PYTHIA [8] does not reproduce the ratio of antiproton to 
proton production seen in the data at any rapidity, while the additional hy- 
pothesis of a baryon junction within the HIJING/B [9] model yields a good 
agreement with both the magnitude and rapidity dependence of the observed 
p/p ratio. 



2 The analysis 

The data presented in this Letter were collected with the BRAHMS detector 
system during 2001. BRAHMS consists of two movable magnetic spectrom- 
eters and a suite of detectors designed to measure global multiplicity and 
forward neutrons [10]. In addition, eight rings of plastic scintillator tiles were 
used to find the collision point and provide a minimum bias trigger [11]. To 
reduce the contribution of background events valid hits in the outer three rings 
were required as part of the offline analysis. Using a GEANT simulation with 
the HIJING event generator [12] as input, it was estimated that this trigger 
setup saw 71±5% of the 41 mb p+p total inelastic cross-section. Spectrometer 
triggers that required hits in several hodoscopes were used in each of the two 
spectrometers to enhance the event sample oi p + p collisions with tracks. For 
this analysis data taken at nine angle settings with respect to the beam were 
used, ranging from 90° to 3° and yielding a rapidity coverage of < y < 3.4 
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for pions. 



Identification of charged hadrons (tt, K, and p) was done primarily through 
time-of-flight measurements. Tracks having a measured inverse velocity ( 
within a ±2cx band of the theoretical value for the appropriate momentum and 
mass, were selected for analysis. In the forward spectrometer where particles 
in general have higher momenta, identification was also provided through the 
recorded radius in a Ring Imaging Cherenkov detector, and via momentum 
dependent cuts in the response of a threshold Cherenkov detector. The de- 
tails of the particle identification and analysis methods used are similar to 
those described in [13,14], but because of the lower particle yield our cal- 
ibrated time-of-flight resolutions are worse than for Au+Au. This mainly 
affects the midrapidity spectrometer, which only has time-of-flight systems. 
For the present analysis a separation of p/K up to p =2.6 GeV/c and K/ir up 
to p =1.2 GeV/c was achieved. 

Charged particle ratios were measured by dividing transverse momentum spec- 
tra, normalized to the minimum bias trigger. By measuring positive and nega- 
tive particles at the same angular setting but with opposite magnet polarities, 
most corrections for geometrical acceptance and detector efficiencies cancel 
out. Figure 1 shows the resulting like-particle ratios as a function of px at 
the extreme measured rapidities of y ~ and y ~ 3. Within our statistical 
errors there is no significant dependence on p^. The ratios were therefore fit- 
ted to a constant over a px range matching the limits of our acceptance (see 
Fig. 1). For most settings this range was 0.5 < pi < 1.5GeV/c, varying by 
< ±0.5GeV/c for the different spectrometer angles. 

The ratios have been corrected for particle absorption and in-flight decay 
as discussed in Ref. [13]. In addition corrections were applied for antiproton 
absorption in the spectrometer trigger slats, which removed ~ 10% of the p 
yield at p <1 GeV/c, dropping to ~ 5% at p =2 GeV/c. Primary particles were 
selected by by requiring the tracks to point back to the beam line, with an 
achieved resolution of a ~ 0.7 cm. For n~ /n + and K~/K + a 3a cut was used, 
while for p/p a 2a cut was set to further eliminate knock-out protons from the 
beampipe. Since the spectrometers have a small solid angle the effects of feed- 
down from weak decays are not large and tend to cancel in the ratios [14] . The 
p/p ratio is exceptional since it is sensitive to the evolution with rapidity of 
the K/p ratio. To estimate the upper limits of this effect, a GEANT simulation 
with published STAR data from p + p collisions y = [15] as input has been 
used. Taking A/p ~ 0.5, assuming a constant behavior with rapidity and that 
A/A ~ p/p ■ K + /K" (see e.g. [16]), the feed-down from A and A were found 
to cause a net increase of p/p at all rapidities. At midrapidity the possible 
contribution is <5%, and at forward rapidity <10%, within our acceptance. 
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Fig. 1. Particle ratios vs. pi at y = (solid circles) and y ~ 3 (open circles). The 
lines show the result of fitting a constant to the data, over the indicated range. The 
shaded area shows our estimate of the systematic error. 



3 Particle ratios vs. rapidity 



Figure 2 shows the resulting ratios of antiparticle-to-particle yields as a func- 
tion of rapidity (left panel). Two independent analyses were performed. By 
comparing these, and by varying both the rapidity and pt intervals, and the 
cuts on the particle identification and projection to the interaction point, our 
point-to-point systematic errors are estimated to be <2% for pions and pro- 
tons, and <3% for kaons. Ratios from measurements with different magnet 
polarities allow us to investigate systematic effects from geometry and nor- 
malization. The combined residual systematic uncertainties from these effects 
and from the absorption corrections are found to be <5%. 

For all three ratios in Fig. 2 there is a clear midrapidity plateau and subse- 
quent decrease with rapidity. The midrapidity values of the ratios are n~ /ti + 
= 1.02±0.01±0.07, K-/K+ = 0.97±0.05±0.07 and p/p = 0.78±0.03±0.06, 
consistent within statistical errors with values extracted from identified par- 
ticle spectra reported by STAR [17]. Numbers at other rapidities are given 
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Fig. 2. Left: Charged particle ratios from p + p at yfs = 200 GeV (solid points) 
compared with Au+Au [13] (open points), and predictions from PYTHIA [8] (solid 
histogram) and HIJING/B [9] (thick dashed line). Right: Ratios shifted by yj,, 
compared with data from NA27 (triangles) at y/s = 27.5 GeV [19]. 

in Table 1. At midrapidity, proton and antiproton production from quark- 
antiquark pairs can be assumed to be identical. Proton excess, defined as 
(N p — Np) / (N p + Np), is therefore due to the transport of baryon number from 
the initial beam. Our p/p ratio would in this interpretation imply a proton 
excess of 12% at midrapidity, carrying baryon number that has been trans- 
ported from the beam region at y — 5.3 [4]. We note that it has been shown 
(see [18]) that one may need to correct for isospin effects before generalizing 
these results from p + p to hadron-hadron collisions, due to the presence of 
neutrons. 

At y < 1.5 the Au+Au ratios for the 20% most central collisions reported 
in [13] are noticeably similar to the present results. Above y = 1.5 the pion 
ratios in p + p start to drop below those for Au+Au and consequently below 
unity, while the kaon and proton ratios remain consistent with the Au+Au 
results over our entire acceptance range. This is surprising in view of the 
different dynamics one might expect for the two systems. A heavy ion system 
has multiple initial collisions as well as significant rescattering and may reach 
thermal equilibrium before freezeout occurs, while the significantly smaller 
p+p system should not interact much beyond the initial reactions. For all three 
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Table 1 



Numerical values for charged particle ratios as a function of rapidity. Errors are sta- 
tistical only. In addition a combined systematic error of 7% for 7r~ /tt + and K~ /K + , 
and 8% for p/p is estimated. 



species the ratios start to decrease above y = 1.5, indicating a transition from 
the string breaking dominated regime at midrapidity to the fragmentation 
region. The drop in the pion ratio at high rapidity can be attributed to isospin 
and charge conservation in the fragmentation region, an effect not seen for 
Au+Au where the high pion multiplicity drives the system towards isospin 
equilibration. 



The right panel of Fig. 2 shows the present data and data from NA27 at s/s 
= 27.5 GeV [19] (open triangles) shifted by the respective beam rapidities. 
Overlaying the two datasets the ratios appear to be independent of the inci- 
dent beam energy when viewed from the rest frame of one of the protons, in 
the region where our rapidity coverage overlaps with that of NA27. This is 
consistent with the idea of limiting fragmentation that has also been observed 
for charged hadrons in nucleus-nucleus collisions [20]. This hypothesis states 
that the excitation of the leading protons saturates at a moderate energy, 
leaving more available kinetic energy for particle production below the beam 
rapidity. We also note a transition in behavior at y — y^ ~ —4, indicative of 
a boundary between the midrapidity and fragmentation regions. Below this, 
at RHIC energies we observe a region of constant relative particle production 
that was not present at \fs = 27.5 GeV. 
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4 Predictions from models 



To interpret these results further, predictions from theoretical models of hadron- 
hadron collisions are confronted with the data. The curves in the left panel 
of Figure 2 compare our results to the predictions of two such calculations, 
PYTHIA Ver. 6.303 [8] 2 and HIJING/B [9], using the same p T range as the 
present analysis. Both models give a good description of the pion data and 
for kaons at midrapidity, but do not reproduce the magnitude of the decrease 
with rapidity seen for K~ / K + as the rapidity approaches that of the frag- 
mentation region. Also, PYTHIA clearly overestimates the p/p ratios. This 
is a well-known problem since PYTHIA employs only quark-diquark break- 
ing of the initial protons, while several authors have pointed out [4,5] that to 
describe stopping at midrapidity in high energy hadronic collisions one needs 
an additional mechanism to transport baryon number away from the beam 
rapidities. 

Based on p + p data from the ISR it has been proposed that other mechanisms 
than quark-diquark breaking, e.g. destruction of the diquark, can transport 
baryon number over a large rapidity range [7]. Subsequently a description 
was formulated of the baryon transport process as arising from gluonic de- 
grees of freedom, with an additional transport component slowly changing 
with incident energy [4]. This can lead to a significant net baryon content 
at midrapidity. Also, data from HERA [21] show a baryon asymmetry, de- 
fined in lepto-production as 2 • (N p — N p )/(N p + N p ), that is significantly 
different from zero. This indicates that baryon transport over 7 units of ra- 
pidity is indeed possible. Together, these theories and observations form the 
basis for implementing the baryon junction [4,22]. This mechanism allows for 
easy transport of baryon number toward midrapidity, while energy balance is 
maintained through an increased production of forward mesons. The baryon 
junction scenario, incorporated as a model prediction in the HIJING/B event 
generator [9] , has successfully predicted the slow y^s dependence of the p + p 
and p + p cross-sections [4]. In Fig. 2 the dashed lines showing the HIJING/B 
prediction for p/p at \/s = 200 GeV, exhibit a much better agreement with 
the data than PYTHIA, both in terms of overall magnitude and the width of 
the distribution. 

In Ref. [23] a baryon junction extension to a quark-diquark breaking model 
of particle production is suggested. It is shown that it is possible to describe 
baryon stopping in p + p and Au+Au collisions using the same parameters for 
the baryon junction couplings, but with different parameter values for SPS 



2 PYTHIA version 6.3 is at the time of writing still labeled as 'experimental', but 
we find no difference in the results between this version and the latest in the 6.2 
series. 
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Fig. 3. y/s dependence of particle ratios at y = (closed symbols) and y ~ 1 (open 
symbols). Circles are the present data, errors are statistical only. Also shown are 
p + p data from ISR (squares) and NA27 (triangles) [19,24]. Solid lines: PYTHIA 
prediction for p + p at y = 0. Dashed lines: same for y = 1. Dotted line in bottom 
panel: HIJING/B prediction for p/p at j/ = 0. 

and RHIC energies. For RHIC, this leads to a prediction that the shapes of 
the rapidity distributions for p + p and Au+Au will be similar for \y\ < 2. The 
similarity shown here of p/p in p + p and Au+Au up to \y\ < 3 supports this 
prediction. 



5 Particle ratio excitation functions 

The present data allow for an extended study of the excitation function of the 
particle ratios around midrapidity. In Figure 3 the present data at y = and 
y ~ 1 are shown, together with fits to ISR data [24] from p+p collisions in the 
range 23 < ^/s < 63 GeV. Where possible the fits have been made over the 
same px range as our data, the notable exception being the p/p ratios at y — 1 
where the ISR data cover 2.0 < p-r <4.0GeV/c. Points from NA27 at y/s = 
27.5 GeV are also shown. Both at midrapidity and at y — 1 the ratios depend 
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Fig. 4. Correlation between K~ /K + and p/p at different rapidities from the present 
data and data at lower energies. The lines show grand canonical model calculations 
for the limit of vanishing strangeness chemical potential fi s = (dashed) and for a 
constant temperature of 170 MeV with unit strangeness saturation [28] (solid). 

logarithmically on y/s, but the slope of this dependence is steeper at y — 1. 
At lower energies there is a significantly larger fraction of K~ and antiprotons 
at y = than at y — 1, but this effect is much smaller at RHIC energies. This 
again indicates that at RHIC there is a midrapidity source that is almost free 
of net strangeness and baryon number. 

The solid and dashed lines in Fig. 3 show the prediction for the particle ra- 
tio excitation function from PYTHIA at y = and y — 1 respectively. At 
midrapidity the ratios are well reproduced at all values of y/s, except for the 
p/p ratio at RHIC energies, but at y — 1 the K~/K + and p/p do not seem 
well described at lower energies. The dotted line shows the prediction for p/p 
from HIJING/B at y — 0, reproducing the result at y/s = 200 GeV but un- 
derpredicting the results at lower energies. For pions and kaons HIJING/B 
reproduces the PYTHIA curves shown. 



6 Ratio correlations over three units of rapidity 



For nucleus-nucleus collisions at ultrarelativistic energies it has been observed 
that almost all particle production ratios can be reproduced by a grand canon- 
ical model description of the emitting source, i.e. with temperature T and 
baryochemical potential fi q as independent parameters [25] . The strange quark 
chemical potential ji s is fixed by conservation of strangeness [26]. In such an 
approach antiparticle-to-particle ratios are controlled by the light and strange 
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quark fugacities, fJ> q /T and fi s /T, respectively, predicting e.g. 

K-/K+ = e 2 ^' T ■ e- 2 ^ T = e 2 ^' T • (p/p) 1 * 



(1) 



For an ideal quark gluon plasma one can expect /jl s — 0, a condition that is 
difficult to achieve for a hadron gas [27]. The analysis in Ref. [13] on data from 
Au+Au collisions at y/s = 200 GeV showed that one can parametrize the kaon 
and proton ratios at different rapidities as a power law: K~ / K + = (p/p) a , with 
a Au+Au _ q 24 _|_ 0.02. Expressing this in terms of chemical potentials gives 
H s ~ 0.28/ig for Au+Au collisions. 

Figure 4 shows a similar analysis based on the present data, where the K~ / K + 
ratios have been interpolated to the same rapidities as the p/p data. A power 
law fit to the present points gives an exponent of a p+p = 0.32 ± 0.04, with 
X 2 /NDF = 1.22. Figure 4 also shows the corresponding results for p + p 
collisions at ^Js = 27.5 GeV at rapidities < y < 3.5, and midrapidity data 
at ISR energies [19,24]. The ISR results are consistent with the power law fit 
to our data, while the = 27.5 GeV data seem to follow a different trend. 

The solid line in Fig. 4 is the prediction of a grand canonical calculation for a 
constant temperature of 170 Me V [28]. This curve gives a good description of 
our Au+Au data, as well as lower energy heavy ion results. For y < 2.0 the 
p + p data are also consistent with this curve, but at more forward rapidities 
they fall below it. Ideally for p + p collisions one would use a microcanonical 
approach in order to exactly conserve quantum numbers in each event. Such 
a description is being developed e.g. by the authors of Ref. [29,30], but they 
also show that the K~ / K + and p/p ratios change by < 4% when going from 
the canonical to the microcanonical description. 

The limit of a canonical ensemble can be reached from a grand canonical 
description by letting all chemical potentials approach 0. In e + + e~ collisions 
such a canonical approach has been successful in describing particle ratios [30], 
but this does not imply that such collisions constitute an ideal quark-gluon 
plasma. Rather it may reflect properties of the hadronization process. In the 
above grand canonical approach, a power law exponent of a = 0.33 implies 
that n s = (see the dashed line in Figure 4 and Equation 1). The fit made to 
the present data suggest that this is the case for all covered rapidities in p + p 
collisions at \/s = 200 GeV. 



7 Conclusions 

In conclusion, the BRAHMS experiment has measured ratios of charged anti- 
hadron to hadron production from p + p collisions at y/s = 200 GeV. All ratios 
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are independent of transverse momentum within errors for pp < 2.0GeV/c. 
For kaons and protons we find an overall consistency with results from Au+Au 
collisions at the same energy over three units of rapidity. The n~ /n + ratio falls 
steadily below the Au+Au results for y = 2.0 — 3.4, as expected from conser- 
vation of initial charge and isospin. When viewed from the rest frame of one 
of the protons all ratios seem to be independent of the projectile beam energy 
over a range of at least one unit of rapidity. Models based on quark-diquark 
breaking of the initial protons give a reasonable description of 7t~/tt + , but 
cannot describe our p/p ratios unless additional mechanisms of baryon trans- 
port are invoked. Introducing a baryon junction scheme to provide additional 
baryon transport to midrapidities yields a good description of our p/p data 
over our full coverage of < y < 2.9. 

After submission we have learned about a midrapidity analysis similar to the 
one presented here, made by the PHOBOS experiment [31]. Their result for 
p/p at y = is somewhat higher than ours, but within errors the ratios 
reported by PHOBOS are consistent with the ones presented in this paper. 
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of Norway, the Jagiellonian University Grants, the Korea Research Foundation 
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